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ABSTRACT 



Context. R Aqr is one of the closest symbiotic binaries and the only D-type system with radial velocity data suitable for orbital 
parameter estimation. 

Aims. The aims of our study are to derive a reliable spectroscopic orbit of the Mira component, and to establish connections between 
the orbital motion and other phenomena exhibited by R Aqr 

Methods. We reanalyze and revise the velocity data compiled by Mcintosh & Rustan complemented by additional velocities. 
Results. We find an eccentric orbit (e = 0.25) with a period 43.6 yr This solution is in agreement with a resolved VLA observation of 
this system. We demonstrate that the last increase in extinction towards the Mira variable in 1974-1981 occurred during its superior, 
spectroscopic conjunction, and can be due to obscuration by a neutral material in the accreting stream. We also show that jet ejection 
is not connected with the orbital position. 
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1. Introduction 

R Aqr is a symbiotic binary surrounded by an hourglass neb- 
ula (Solf & Ulrich [19851 ). and both its cool and hot components 
are unique. The cool component is a Mira variable with a pulsa- 
tion period of 3ST^ with SiO and H2O maser emission, which is 
rare, with only three among 48 symbiotics with Miras exhibit- 
ing this emission (Whitelock 120031 ) . The hot component shows 
sporadic, jet ejections (e.g. Kellogg et al. |2007| and Nichols et 
al. 12007b . and Nichols et al. (12007b also detected a 1734 s peri- 
odic oscillation in Chandra X-ray observations, which suggested 
that the hot companion is a magnetic white dwarf. Willson et al. 
( 11981b interpreted a brightness reduction in pulsation maxima 
as the eclipsing of the Mira by a companion surrounded by an 
extended gas cloud. This interpretation was supported by two 
facts. First, during the last hypothetical eclipse (1974-1980), a 
minimum was observed in the near-IR light curve similar to that 
corresponding to dust obscuration in other symbiotic Miras (e.g., 
Whitelock et al. |1983l ). Secondly, a spectrum obtained during the 
first such event (1886-1894) exhibits only faint emission lines of 
hydrogen and no trace of M type spectral features (Townley et 
al. [1928). In the historical, visual, light curve, three phases of 
reduced maxima are present spaced at 44 yr intervals and this 
value was adopted by Willson et al. ( 11981b as the orbital period. 
Additionally, the O - C diagram showed variation with a period 
of ~22 yr. Wallerstein ( 1986) used a number of spectral absorp- 
tion lines associated with the Mira in an attempt to determine the 
orbital elements of R Aqr. His results were consistent with the 
44 yr period, but uncertain. Hinkle et al. ( 11989b combined these 
archival data with radial velocities (RV) determined from near- 
IR CO and Ti I lines, and obtained a 44 yr period and a highly 
eccentric, e = 0.6, orbital solution. 



Mcintosh & Rustan (12007b collected radial velocity data of 
SiO maser emission of R Aqr (Lepine et al. 119781 Cohen & 
Ghigo[l98 0l Spe ncer et a l.fTMTl Lane[l982l Hall et al. [19901 
Ch o et a l. 1199^ Bob oltz 119971 Hollis et al. [200 01 Alco lea et 
al. [T999I Pardo et al. 120041 Mcintosh & Rustan 120071 ) com- 
plemented by previous radial velocities derived from optical 
absorption lines (Merrill 119501 Jacobsen & Wallerstein 119751 
Walle rstein 119861 ) and near-IR CO and Ti I lines (Hinkle et al. 
119891) . All of these features are associated with the cool compo- 
nent or material in its close neighborhood and reflect its motion 
around the mass center. In astrophysical masers, emission orig- 
inating in many maser spots moves around the Mira variable. 
Spectra of masers exhibit a series of emission lines with vari- 
ous velocities associated with different maser spots. Mcintosh 
(120061) studied changes in the velocity centroid (VC) of SiO 
maser emission from 76 single evolved giants, and showed that 
the mean difference in the VCs measured at different epochs was 
only 0.065 ± 2.00km s"', and that the VC of the emission dis- 
tribution can be used to derive the velocity of the star. Based on 
this result, Mcintosh & Rustan (12007b adopted the VC in their 
work for all kinds of velocity data analysis. 



Send offprint requests to: M. Gromadzki 



Mcintosh & Rustan (12007b also estimated orbital elements 
using a program developed by Gudehus (120011 ). They inferred 
an eccentric orbit, e - 0.52, with a period 34.6 yr. Their orbital 
period was significantly shorter than any previous estimates as 
well as the period derived in our present study. Their other were 
of reasonable accuracy if we interpret their system mass, equal 
to 0.043 M0, as the mass function, and their semi-major axis 
of the system, of 3.7 AU, as the semi-major axis of the orbit 
of the Mira component. We also note that not all points listed 
in their Table 1 were plotted on their figures showing VC versus 
orbital phase, namely, they did not measurement from Alcolea et 
al. (11999b . In the text, the authors wrote that points from Hollis 
et al. (12000b are listed in their Table 1, although they were not 
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displayed there. These data points were plotted on the figures 
showing VC versus Julian days, but they were missing on the 
figure of VC versus orbital phase. The authors also did not com- 
pare the orbit solution with the resolved image of the binary 
(HoUis et al. |1997| l. nor include radial velocity measurements be- 
fore 1946 (Merrili ri935l|1950] l. These points were also ignored 
in all previous papers about the orbit estimation from radial ve- 
locities (Wallerstein |1986l Hinkle et al. |1989l l. The authors also 
neglected radial velocity measurements of the SiO maser by 
Zuckerman et al. SWM . H ollis e t al. ( fT986l l. H all et a l. ( fT987] l. 
All en et a l. ( fT989] l. He ske ([19891 ). Jewell et al. ([1991] ). Patel et 
al. fmTh Gray et al. ([T995]l, Sch warz et al. ([1997), HolHs et 
al. ((2 000 , 200 Ij Imai et al. (ISOOTT l. Cotton et al. (12004, 20061 
and Kang et al. ( 120061 ). We also note that the observations stud- 
ied by Pardo et al. ( 120041 ) included the observations analyzed by 
Alcolea et al. ( 119991 ). which included the data used by Martinez 
et al. (119881 ). These three data sets are thus not entirely different, 
and the later studies were based simply on data collected over 
longer time. 

This paper contains an analysis of the revised radial velocity 
data listed in Table 1 of Mcintosh & Rustan (2007) and addi- 
tional velocity data (listed here in Tabled]), as well as discussion 
of the connections between the orbital motion and other phe- 
nomena shown by R Aqr 

2. Orbit determination 

The data from Table 1 of Mcintosh & Rustan ( 120071 ) were re- 
vised and supplemented by additional radial velocity measure- 
ments from Merrill ( 119351 11950] ). Zuckerman ( 119791), M artinez 
et al. ([T988] ) Jewe ll et al. ([T99T] ), Schwarz et al. ([79^5] ), HolHs 
et al. ( |2000l [2001] ). Cotton et al. ([2004] [2006]) . and Kang et al. 
([2006] ). 

The blue-violet absorption lines in Mira variables show con- 
stant shift of between -5 km s"' and -10 km s"' with respect 
to near-IR velocities (Hinkle et al. |1984t and references therein). 
In our study, a shift equal to -6 km s"^ was adopted. The blue- 
violet absorption lines do not show any variation with pulsation 
period and do not need any additional correction. Absorption 
lines from the red part of the optical spectrum exhibit signifi- 
cant variation but there is no way of correcting for this effect 
in a single measurement. For this reason, we did not use three 
data points from Jacobsen & Wallerstein ( 119751 ), which had been 
previously used for the spectroscopic orbit determinations by 
Waller stein ( [1986] ), Hinkle et al. ( [TM9] ), and Mcintosh & Rustan 
([2007] ). 

In the near-IR, the radial velocity (RV) variations of R Aqr 
were dominated by stellar pulsations (AV ~ 40 km s"'). To re- 
move changes due to the pulsations, Hinkle et al. ( 119891 ) com- 
plete a linear fit to the infrared velocity curve between phase 
0.1 and 0.8, and then adopted as the stellar velocities the resid- 
uals remaining after subtracting the velocity-curve fit from each 
measurement plus the center-of-mass velocity of the Mira of 
-29.1 kms"'. The value of the center-of-mass velocity was es- 
timated from the CO Av = 3 velocity at phase 0.36, following 
the results of Hinkle, Scharlach & Hall ( I1984D . Unfortunately, 
Hinkle et al. (19891) used an old ephemeris of Campbell (1955]). 
Since R Aqr displays significant variation in its pulsation period 
(see Fig. [1]), we repeated these corrections with the ephemeris 
JD(max) = 2 442 404.2(+0.4) -H 388. 1 (±0 1 ) x £ , derived for the 
1975-1989 period covering Hinkle et al. ( 119891 ) observations. 

To calculate the orbital solution all data points were averaged 
in 387-day bins corresponding to the pulsation period of the Mira 
component. The velocity changes due to the orbital motion are 



negligible during such an interval. The data obtained from differ- 
ent methods were binned separately. Since different maser tran- 
sitions are formed at different distances from the stellar surface, 
they were also binned separately. We excluded single measure- 
ments for various maser transitions (Allen et al. 119891 Patel et 
al. [T992] Gray et al. [T995] Imai et al. HQOI) obtained for epochs 
covered by more extensive data sets. In Tabled we list the radial 
velocity data used to derive the orbital elements of R Aqr and 
then converted to local standard of rest (LSR) velocities, which 
are also plotted in Fig. [1] 

We also decided to exclude from our analysis the SiO maser 
data obtaine d duri ng the period JD2446 000-47500 (data from 
Holliset al. [T986] Heske [T989] Alcolea et al. [T999] and Hall et 
al. 119901 ) because the maser emission showed at that time signif- 
icantly different behaviour e.g. a few fainter emission, irregular 
contours than in later epochs, as well as the masers in isolated 
Miras (see Pardo et al. 120041 ). The SiO masers in R Aqr display a 
ring-like morphology ~ 31 mas (~ 6.2[^f/200pc] AU) in diam- 
eter (Boboltz et al. 119981 ). So, these masers lie very close to the 
stellar surface, namely at a distance of only ~ 1 .5 radii from the 
Mira (see Sect. 3) in a region dominated by stellar pulsations and 
permeated by circumstellar shocks. In addition, the SiO masers 
can be seriously affected by tidal interactions during the perias- 
tron passage when the size of the SiO maser region became com- 
parable to the Roche lobe radius of the Mira. In fact, according 
to our orbital solution (see below) the Alcolea et al. ( 119991 ) ob- 
servations were obtained around periastron. The deep minimum 
in both VCs and flux of maser emission that appeared around 
JD2 446 700 ^7 300, coincided with the jet formation (Nichols 
et al. 120071 ) and a rise in the UV continuum (Meier & Kafatos 
119951 ). We suspect that the SiO masers may have been partly dis- 
rupted by rising UV radiation during the jet formation, and tidal 
interaction during the periastron passage and then gradually re- 
build. 

To identify the spectroscopic orbit, we first estimated a pe- 
riod of RV variations of 43.6 yr using Ortfit (based on the method 
of Schwarzenberg-Czerny 1996), and then calculated the orbital 
parameters: Tq, the time of periastron passage, y, the baricentral 
velocity, K, the semi-amplitude, w, the longitude of periastron, 
and e, the eccentricity with fixed period. We applied Bertiau's 
program ( 119671 ) based on the Lehmenn-Filhes method. As a re- 
sult, we obtained an eccentric, e - 0.25 ±0.07, orbit with reason- 
able errors (Table|2]). Figure [2] displays the data and the best-fit 
velocity curve versus orbital phase. 

Figure [1] also compares the RV variations with the O - C 
data for the Mira pulsations based visual data collected by am- 
ateur observers from AAVSO. The 0-C is a plot showing the 
difference between the observed (O) times of maximum light 
and the values calculated according to an adopted ephemeris (C) 
as a function of time. A description of the 0-C method can be 
found in a review by Zhou (120031 ). To derive a new pulsation 
ephemeris, we first estimated the time of each maximum in the 
light curve by fitting a third-order Fourier polynomial. Then a 
linear ephemeris: 

JD(max) = 2 416 070 + 4 + 387.30 + 0.07 x E 

was fitted to the obtained times of maxima. The resulting 0-C 
diagram shows a parabolic shape that corresponds to an increase 
in the period rate of 0.0 1 3 ±0.002 day per cycle. After subtracting 
this parabolic trend, a sinusoidal variation with a period of 22.5 
years and amplitude of about 8 days (Fig.[T]) is visible. 
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Fig. 1. The O - C diagram (top) and RV variations (bottom) for R Aqr. Observations from different sources are marked like in 
Table[T] The solid curve represents the fitted orbital elements (Table[2]i. The dotted curve represents the orbital model proposed by 
Hollis et al. ( |1997| i. The dashed curve represents solution, with modified w = 87 which better fits the VLA observations. Bars show 
times of hypothetic ecUpse, short arrows represent moments of jet ejection (Nichols et al. l2007l l. whilst long arrows show periastron 
passage. 



Table 2. Orbital elements for R Aqr. 



Element 


Value 


Porb (day) 


15 943+471 


y (km s"') 


-24.9±0.2 


ifCkms"') 


4.0±0.4 


e 


0.25+0.07 


CD (deg) 


106+19 


To 


2444 019+728 


S(0 - cf 


40 


a (km s"') 


1.11 



3. Discussion 

Our orbital solution yields the semi-major axis of the Mira of 
KP 



flg sin i = — Vl -e2 = 5.68![J;g AU, 
and the mass function of 

/(M) = — ■ PK^d - e^)^/^ = ^^"'^"'^^ = 0.096 t^t?Mo, 
■'^ 27rG (Mh-HMg)2 ° 

where Mg, and Mj, represent masses of the giant and the hot 
component, respectively, and the errors are set by the uncertain- 



ties in the orbital parameters (see Table |2]i. The mass function 
sets a lower limit to the hot companion mass. 

For further discussion, we assume that the symbiotic binary 
R Aqr consist a Mira variable and a white dwarf (WD) compan- 
ion with an orbit inclined at ; = 70° (~ 72° according to Solf & 
Urlich [T985] l. Assuming plausible components masses, Mg = 1- 
1 .5 Mo, and Mh =0.5-1 .4 Mq, for the Mira and the WD, respec- 
tively, the total mass of system M,ot = Mh + Mg should be within 
the range 1.5-3 M©. The component masses are also limited by 
the orbital solution. Therefore, if we consider /(M) - 0.064- 
0.138 Mo and Mg = 1-1.5 Mq, we derive Mto, = 1.6-2.5 Mo, 
Mh = 0.57-1.02 Mo, and the mass ratio q = Mg/Mh = 1.2-2.1 
(Fig- ID- These values are in good agreement with the commonly 
adopted model of symbiotic binaries. The semi-major axis of the 
system, a, is 14.2-16.8 AU, which corresponds to 71-84 mas on 
the sky for a distance of d - 200 pc (see below). 

Our orbital solution can be tested and refined by using re- 
solved observations of R Aqr (Hege et al. 119911 Hollis et al. 
[1997)- In particular, the binary components of R Aqr were re- 
solved by VLA observations (Hollis et al. |1997t obtained on 
20 Nov 1996 (JD=2 450407). The measured separation between 
the VLA image in SiO maser v=l, 7=1-0 transition and the 
continuum emission in a 50 MHz bandwidth at 43.165 GHz as- 
sociated with an Hii region surrounding the hot companion was 
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Table 1. Radial velocities of R Aqr used for orbital elements determination. 



JD 


RV [kms-'l 


Spectral rang 


e 


Symbol 


References 




2422255 


-21.6 


Visual 




C 


Merrill 


1935 




2423039 


-22.0 


Visual 




C 


Merrill 


1935 




2429581 


-30.7 


Visual 




a 


Merrill 


1950 




2430336 


-30.2 


Visual 




a 


Merrill 


1950 




2432235.5 


-28.5 


Visual 




a 


Merrill 


1950 




2433019.3 


-29.0 


Visual 




a 


Merrill 


1950 




2441237 


-21.7 


Visual 




• 


Jacobsen & WallersteinI 19751 


2442939 


-22.0 


Near-IR 






Hinkle etal.|1989l 


2443527.5 


-26.6 


V = l,y = 1 - 





+ 


Lepine et al.ll9781 ZuckermanI 19791 


2444078 


-25.9 


V = i,y = 1 - 





■ 


Cohen & Ghieol 19801 Spencer et al.|1981| Lanel 19821 


2444356 


-27.0 


V = i,y = 1 - 





X 


Lane|l982| 




2444509 


-25.6 


Visual 




o 


WallersteinI 19861 


2444748 


-30.1 


Visual 




o 


WallersteinI 19861 


2445535 


-27.7 


Near-IR 






Hinkle etal.ll989l 


2445574.5 


-27.6 


V = i,y = 1 - 





▲ 


Cho et al.|1996l Jewell et al.|1991| 


2445862 


-30.1 


Visual 




o 


WallersteinI 19861 


2445890.6 


-28.0 


Near-IR 






Hinkle etal.[l989| 


2446378.5 


-27.8 


Near-IR 




5« 


Hinkle etal.n:989| 


2447247 


-28.5 


V = i,y = 1 - 





A 


Martinez etal. 1988] 


2447338 


-28.9 


Near-IR 




^ 


Hinkle et al. 1989) 


2447634 


-28.5 


V = i,y = 1 - 





A 


Martinez etal.fMS] 


2447870 


-27.1 


V = i,y = 1 - 





T 


Pardo et al.'2004| 


2448240 


-27.0 


v= l,y = 2- 


1 


♦ 


Schwarz et al. 19§5\ 


2448283 


-28.3 


V = i,y = 1 - 





T 


Pardo et al. 2004) 


2448696 


-27.4 


v= i,y= 1- 





T 


Pardo et al. 2003 


2448998.4 


-26.4 


V = i,y = 2- 


1 


♦ 


Schwarz et al. 19pg] 


2449435 


-27.0 


V = i,y = 1 - 





T 


Pardo et al. 2004| 


2449439 


-26.3 


v= i,y = 2- 


1 


♦ 


Schwarz et al. I9S] 


2449804 


-27.0 


V = i,y = 1 - 





T 


Pardo et al. 2004| 


2450068.8 


-26.8 


V = 1,7 = 1 - 








Boboltz et al. 19^1 


2450407 


-24.0 


V = i,y = 1 - 





□ 


HoUis et al. 12000] 


2450948 


-24.0 


V = i,y = 1 - 





□ 


Hollis et al.l2000l 


2451390.4 


-25.6 


v= 1,7 = 2- 


1 


o 


Kana et al. 120061 Hollis et al. 120001 


2451785.8 


-23.2 


v= 1,7 = 2- 


1 




Kana etal. 120061 


2451892 


-24.3 


v = i,y = 1 - 





e 


Hollis et al. 120011 Cotton et al.l2004l Mcintosh & Rustanl2007l 


2453253 


-23.0 


V = 1,7 = 1 - 





Kl 


Cotton etal. 120061 


2453784.3 


-22.7 


V = 1,7 = 1 - 





V 


Mcintosh & Rustan 2007] 


2454112 


-22.8 


V = 1,7 = 1 - 





V 


Mcintosh & Rustan i2007l 



p - 55 + 2 mas (ll[c//200pc] AU), and the position angle was 
6 = 18° + 2°. This can be compared with the projected com- 
ponent separation on the sky plane calculated from our orbital 
elements. The orbital elements (Table |2|i combined with the or- 
bital inclination of / = 70° (Solf & Urlich [T985] l inferred the 
projected separation of p - 0.44a ~ 31-37 [c//200pc]"' mas, 
and the difference between the apparent binary position angle 
and the position angle of the line of nodes (the binary orienta- 
tion on sky) of - Q. = -109°. 

Our predicted component separation on the sky is consis- 
tent with the VLA observations only if the distance is somewhat 
lower than 200 pc or we underestimate the p value due to uncer- 
tainties in our orbital elements. 

The predicted value of p should in fact increase, changing co 
(see also Fig. |4]i, and/or increasing e, and/or decreasing Tq. For 
example, p should increase to 0.65a (46-55 [^f/200pc]"' mas) 
just by adopting a lower value of to - 87, i.e. lower only by 
an amount approximately equal to the standard-deviation uncer- 
tainty of O). 

The distance to R Aqr is also relatively uncertain, and pub- 
lished estimates based on various methods are in the range 180- 
26 pc ( e.g. 180 pc by Solf & Ulrich 1985, 181 pcby Lepine et 
al. |1978l 240 pc using the most recent period-luminosity relation 
from Whitelock et al. |2008l and 260 pc by Baade [T943] [T944l l. 



The Hipparcos parallax of R Aqr is 5.07±3.15 mas (Ferryman 
et al. 119931 ). which corresponds to a distance of 197^^^^^ pc. Our 
determination of the parameters of the spectroscopic orbit is thus 
in a good agreement (within < l-cr errors in the orbital parame- 
ters) with those inferred from the VLA resolved observations. 
The value of 6* = 18°, measured with the VLA, implies that 
Q. = 127°. Thus, the position of the binary orbit on the sky repro- 
duces well the general picture of the nebula around R Aqr, i.e. 
the orbital plane is perpendicular to the jets (see Fig. 7 in Solf & 
Urlich |T983] l. 

The orbital parameters derived for the speckle observations 
(Hege et al. |1991| l completed on 16 Oct 1983 are inconsistent 
with our orbit in any case, because the measured separation, p = 
124 ± 2 mas which coiTesponding to 25[(i/200pc] AU, is too 
large for any orbital phase. Hege et al. (|1991| l probably did not 
detect the Mira companion but rather an He emission region in 
the SW counterjet. 

Our orbital solution predicts the times of spectroscopic con- 
junctions at Tconji = 2450 859 (Feb 1998), reonjii = 2443 658 
(May 1978) for the inferior ( the Mira is in front of the hot 
component), and the superior (the Mira behind) conjunction, re- 
spectively. The superior conjuction coincides with the last dust 
obscuration/eclipse (see also Fig. |4]i. Although the conjunction 
times could differ by up to 1-2 years due to uncertainties in the 
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Fig. 2. Bottom: The RV data and fitted orbital curve (Table |2]i 
against orbital phase. Symbols are the same as in Fig[T] The solid 
curve represent our orbital solution. The dashed curve represents 
solution, with modified w = 87 which better fits the VLA obser- 
vations. Top: Difference between the RV data and our orbital fit 
(solid line in the lower pannel). 
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Fig. 3. The permitted component masses constrained by the ob- 
served mass function as well as its maximum and minimum val- 
ues (solid lines labeled with the values of f{M)), and assuming 
i = 70°. Dotted and dashed lines show the solution for / - 75°, 
and i = 90°, respectively. 

orbital parameters, especially Tq, e, and oj, the obscuration lasted 
about 6 years, and the coincidence between these two events is 
obvious. 

Adopting plausible binary parameters, i.e. the mass ratio 
q ~ 1.6, Mtot = 2 Mq, and a = 15 AU (see above), we esti- 
mate the minimum component separation to be a(l -e) ~ 1 1 AU, 
and the Roche lobe radius for the Mira, to be /?l ~ 4 AU. The 
near-IR interferometry of R Aqr gave the diameter of the Mira 
expressed as a diameter of uniform disk, @uo = 14.06-20.8 mas 
in K, 17.7 mas in /, and 17.7-19.07 mas in H, in intermedi- 



ate and minimum pulsation phase, respectively (van Belle et al. 
[T996t Millan-Gabet et al. [20551 Ragland et al. |2506l ). These val- 
ues correspond to the average radius of the Mira component, 
Rg ~ 2 AU. Even during periastron passage, the Mira variable 
therefore remains relatively far from filling the Roche lobe. 

Many detached interacting binaries, including the symbiotic 
binaries as well as many X-ray binaries show evidence for much 
higher mass transfer rate than predicted by models of spherically 
symmetric winds. This phenomenon can be accounted for by 
wind focused towards the compact component. Podsiadlowski 
& Mohamed ( I2007I I proposed a model for o Cet binary sys- 
tem that can explain this focusing. In their model, a slow wind 
from Mira fills its Roche lobe and then the matter falls - via the 
Li Lagrangian point - into an accreting stream onto an accre- 
tion disk around the companion. This model is even more rele- 
vant to R Aqr where the component separation is significantly 
smaller than that in o Cet. In particular, it could explain the in- 
crease in extinction towards the Mira during the last "eclipse" in 
1974-1981 caused by a neutral material in the accreting stream. 
According to our orbital solution, the streaming material can ob- 
scure the Mira at that time, and in addition, the proximity of 
periastron (see Fig.|4| could give rise to enhanced mass transfer. 
A similar enhanced, wind-obscuration scenario for R Aqr was 
discussed by Mikolajewska & Kenyon ( |I9921 l, although in that 
case the source of the increased mass loss was proposed to be 
a helium flash above a massive core (Mcore ~ I.36M0) of the 
Mira. 

The orbital configuration of R Aqr could also explain the ob- 
servations of the OH and H2O masers that originate in circum- 
stellar regions more distant from the Mira variable (e.g. H2O 
masers zone is 10-20 AU) than the SiO maser emission region. 
The OH and H2O maser lines were not detected in 1979 and 
1984 (Cohen & Ghigo [19801 Norris et al. [T98l) whe n the hot 
companion was in front of Mira. Ivison et al. ( |I994| l reported 
detection of weak OH and H2O maser emission in 1993, when 
the Mira was in front of the hot companion and the masers were 
shielded from its radiation. They claimed, however, that the OH 
maser detection was only tentative, and that the velocity mea- 
surement was very unusual (an expansion velocity of 10 km s"' 
is rather large for this system), and only the red peak should 
have been detected, and additional spectra obtained using the 
VLA during 1994 June did not show the OH maser line (Ivison 
et al. I1998I I. The H2O line appeared more promising, although 
weak for a Mira variable with relatively strong SiO masers. It 
was also clearly detected during 3 epochs in 1993-1995 (Ivison 
et al. 119981 ). Ivison et al. d 1 9981 ) proposed three possible mech- 
anisms that could inhibit the OH and H2O maser formation in 
their natural environment, involving either the orbital motion of 
the companion or its UV radiation and fast wind, or line obscu- 
ration at low frequencies by optically thick, ionized gas. 

The apparent correlation between the orbital position and the 
0-C for the Mira pulsations is extremely interesting. In particu- 
lar, minima in the (9-C diagram occur around the periastron (see 
Figd)- This may be due to some distortion of the Mira caused 
by tidal interaction during the periastron passage, for example 
if the Mira became elongated towards the hot component, the 
pulsation may need more time for propagation in this direction. 
According to our orbital model, the main maximum of (9 - C 
appears close to the inferior conjunction, when the Mira eclipses 
the hot companion, whereas the secondary, lower maximum ap- 
pears in the middle of the two conjunctions, when both compo- 
nents are well separated and we can see the elongated side. 

Finally, the orbital solution proposed is inconsistent with the 
17 yr orbital period for R Aqr proposed by Nichols et al. ( |2007t . 
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Fig. 4. The orbit of the Mira and hot (x) component in the 
R Aqr binary system in steps of = 0.1. In this representa- 
tion, the stars move anti-clockwise. The dotted circle represents 
the Mira boundary at <p = (periastron passage). The solid dot 
marks the mass center Axes are in units of the semi-major axis 
a. The positions of Mira at the beginning and end of eclipse 
is marked by El, and E2, respectively whereas the main and 
secondary maxima in (9 - C diagram are denoted by Ml and 
M2, respectively. Jl and J2 mark the two jet ejection episodes. 
Open squares connected with dotted line and marked by V show 
the component position during the VLA observation (Hollis et 
al. |1997l l. Open dots connected with dotted line and marked by 
S show the component position during the speckle observation 
(Hege et aL lTWIT i. 

They noticed that radio observations obtained in 1987 and 2004 
exhibited non-thermal spectra, which are usually interpreted to 
be a signature of a newly emitted jet (in stellar sources older jets 
are generally thermal), suggested that the jet formation is cou- 
pled with periastron. In our opinion, such a short period appears 
unrealistic not only because it is not reflected by the RV variation 
but also because in that case there would not be sufficient room 
for a Mira with a dust envelope. The 17-yr period implies a semi- 
major axis of ~ 8 AUa whereas the minimum component sepa- 
ration for D-type symbiotic binaries is about 10-15 AU, which 
for typical masses (see above) indicates orbital periods longer 
than 20 yrs (Mikolajewska |T999] l. It appears that the jet ejection 
in R Aqr connot be connected with the orbital motion. 

4. Conclusions 

Based on published radial velocities of its Mira component, we 
have derived new orbital parameters for the symbiotic binary 
R Aqr. We have found in particular, the orbital period of 43.6 yr, 
and showed that the mass function is consistent with the pres- 
ence of a typical, 1-1.5 Mq, Mira variable accompanied by a 
0.6-1 Mq white dwarf. We also showed that our spectroscopic 
orbit is consistent with the VLA astrometry (Hollis et al. |19971 i 
and inconsistent with the speckle interferometery (Hege et al. 
|1991t . Our orbital model allow us to interpret the "eclipses" of 



Mira as obscuration by the accretion stream. We also showed 
that the jet ejection (Nichols et al. 120071 ) is not connected with 
the orbital position. Finally, we note that September 2012 should 
be a perfect time to complete resolved imaging of R Aqr, since 
the Mira component will then be at minimum brightness and the 
separation between the components on sky will be close to its 
maximum possible value. 
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